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Abstract 

The b quark forward-backward asymmetry has been measured using hadronic decays col- 
lected by the OPAL experiment at LEP. TP — > bb decays were selected using a combination of 
secondary vertex and lepton tags, and the sign of the b quark charge was determined using an 
inclusive tag based on jet, vertex and kaon charges. The results, corrected to the quark level, are: 

= 0.0582 ±0.0153 ±0.0012 at 89.50 GeV 

A^B = 0.0977 ±0.0036 ±0.0018 at 91.26 GeV 

= 0.1221 ±0.0123 ±0.0025 at 92.91 GeV 

where the first error is statistical and the second systematic in each case. Within the framework of 
the Standard Model, the result is interpreted as a measurement of the effective weak mixing angle 
for electrons of sin^ ff'^j''' = 0.23205 ± 0.00068. 
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1 Introduction 

The measurement of Ap^, the forward-backward asymmetry of b quarks produced in e'^e" — > bb 
events, provides an important test of the Standard Model, allowing the effective weak mixing angle 
sin^ 0^ to be determined with high precision The differential cross-section for the production of 
bb pairs can be written as 

da 9 8 k 

1 + cos" 9 + -A^r, cos a , 



d cos 6 3 

where 9 is the angle between the directions of the incoming electron and outgoing b quark, and where 
initial and final state radiation, quark mass and higher order terms have been neglected. The Standard 
Model prediction of the Z*' pole asymmetry A^'^ can be written 



where gy' and (f^ are the effective vector and axial- vector couplings of the electron and b quark to 
the Z*^. The effective weak mixing angle sin^ ^^'^ ^ charged fermion f can be expressed as 




where qi is the electric charge of the fermion in units of the electron charge. With the values of 
the electron and b quark couplings predicted in the Standard Model, the asymmetry is mainly 
sensitive to the weak mixing angle for electrons, sin^ 9^'^ ^ and insensitive to that for b quarks. The 
weak mixing angle sin^ 9'^''^ can therefore be determined from the measured asymmetry within the 
context of the Standard Model, which also predicts the centre-of-mass energy dependence of the 
asymmetry arising from Z'^-7 interference |||]. 

The most sensitive measurements of sin^ ^w'*^ LEP come from the measurements of the b 
quark asymmetry, using techniques based on jet charges, secondary vertices and high momentum 
leptons |,|,|]. This paper reports an improved measurement of using jet, vertex and kaon 
charges combined in an inclusive tag. Compared to the previous OPAL analysis using jet and vertex 
charges Q, the measurement is improved by using the high-performance b-tagging technique developed 
for the measurement of i?b incorporating both vertex- and lepton-based b-tags, b quark charge 
tagging methods developed for B*^ oscillation and CP-violation measurements an increased 

angular acceptance and a more sophisticated fitting technique measuring more of the required event 
properties from the data themselves. The data sample is also increased by adding about 0.5 million 

decays recorded primarily for calibration purposes during the LEP2 physics programme between 
1996 and 2000. 

A brief outline of the analysis method is given in the following section, followed by a description 
of the data sample, the bb event tagging and the b quark charge tagging in sections ^ to ^. The fit 
method and results are described in section ^, and a discussion of systematic uncertainties is given in 



section 
section 18 



7. A summary of the asymmetry results and interpretation in terms of sin ^w''^ are given in 



2 Analysis method 

The analysis method exploits the structure of Z^ bb events, which tend to be composed of two 
back-to-back jets, each containing the decay products of one of the b quarks. Each event was divided 
into two hemispheres by the plane perpendicular to the thrust axis and containing the interaction 
point, and the two hemispheres, each typically containing one b jet, were considered independently. 
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The direction of the thrust axis, in particular its polar angle Ot, was used as an estimate of the 
original b quark direction[|. The hemisphere containing the positive z-axis, i.e. the outgoing electron 
beam direction, was labelled 'forward', and the other hemisphere labelled 'backward'. Two b-tagging 
algorithms, based on secondary vertices and high momentum leptons, were applied to each hemisphere, 
and used to define four classes of b-tagged hemispheres of differing purity (see section ^ . The numbers 
of events with b-tags in one or both hemispheres, together with externally- input values of -Rb and Rc 
(the fraction of hadronic decays to bb and cc) provide enough information to determine the b and 
c quark tagging efficiencies for each tag with only small dependence on Monte Carlo simulation. Jet, 
vertex and kaon charge information in each hemisphere was then used to determine the production 
flavour of the b hadron in the hemisphere, and hence the sign of the underlying quark charge (b or 
b — see section |5|). The two hemisphere determinations were combined to produce a better estimate for 
the whole event. Since each event contains a b and a b quark, the fraction of events where the two tags 
agree can be used to produce an estimate of the reliability (mistag fraction) of the production flavour 
tag. The asymmetry Ap-^ was then extracted from the production flavour tag distributions in forward 
and backward hemispheres, for five bins in |cos^r| and 14 event classes with different combinations 
of b-tags in each hemisphere. 

In this method, the most important quantities needed for the analysis, i.e. the b- and c-tagging 
efficiencies and the fraction of hemispheres with an incorrect production flavour tag, are extracted 
directly from the data as a function of |cos0t|. Monte Carlo simulation is needed to determine 
the tagging efficiencies and mistag fractions for light quark events, the mistag fraction for charm 
events, and the effects of correlations between the two hemispheres in bb events, which result in 
small corrections to both the b efficiencies and mistag fractions. The uncertainties in all these input 
quantities result in systematic errors that are much smaller than the data statistical error. 



3 Data sample and event simulation 

The OPAL detector is well described elsewhere 10, 11]. This analysis relies mainly on charged 
particle track reconstruction using the central tracking chambers and the silicon microvertex detector. 
The latter was first operational in 1991, providing measurements in the r-(p plane only. In 1993 it 



was upgraded to measure tracks in both r-(/) and r-z planes |l^, and in 1996 the cos 6 coverage for at 
least one silicon measurement was extended from | cos9\ < 0.83 to | cos9\ < 0.93 |11|. To account for 
the changing detector performance with time, the analysis was performed separately for data taken 
in 1991-1992, 1993, 1994, 1995 and 1996-2000, and the results were finally combined. 

Hadronic decays were selected using standard criteria, as in The thrust axis direction was 
calculated using charged particle tracks and electromagnetic calorimeter clusters not associated to any 
track. The polar angle of the thrust axis 6t was required to satisfy |cos0t| < 0.95. The complete 
event selection has an efficiency of about 95 % for hadronic decays and selected 3 755 967 data 
events. Of these, around 5 % were recorded at centre-of-mass energies approximately 2 GeV below the 

peak, and 7 % approximately 2 GeV above the peak. This allows the b quark forward-backward 
asymmetry to be measured precisely at three separate energy points. 

Charged tracks and electromagnetic calorimeter clusters with no associated track were combined 



into jets using a cone algorithm |12] with a cone half-angle of 0.65 rad and a minimum jet energy of 
5 GeV. Using a cone rather than a recombination based algorithm increases the fraction of tracks in 
the jet coming from the b hadron decay, which improves both the b-tagging and production flavour 
tagging performance. The transverse and longitudinal momenta of each track were defined relative 
to the axis of the jet containing it, where the jet axis was calculated including the momentum of the 
track. 



right-handed coordinate system is used, with positive z along the e~ beam direction and the a;-axis pointing 
towards the centre of the LEP ring. The polar and azimuthal angles are denoted by 6 and (j>, and the origin is taken to 
be the centre of the detector. 
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Monte Carlo simulated events were generated using JETSET 7.4 [13| with parameters tuned by 



OPAL |14|. The fragmentation function of Peterson et al. |15| was used to describe the fragmentation 
of b and c quarks. The generated events were passed through a program that simulated the response 
of the OPAL detector |16| and through the same reconstruction algorithms as the data. 



4 Tagging bb events 

Two methods were used to tag bb events, based on displaced secondary vertices and high momentum 
leptons. The first method exploits the long lifetime, hard fragmentation, high decay multiplicity 
and high mass of b hadrons, and is fully described in The primary vertex position was first 
reconstructed separately in each event hemisphere, using tracks from that hemisphere combined with 
a common beamspot constraint. Reconstructing separate primary vertices in each hemisphere strongly 
reduces inter-hemisphere tagging correlations. An attempt was then made to reconstruct a secondary 
vertex in each jet of the event, using a subset of well-measured tracks with momentum p > 0.5 GeV. 
If a secondary vertex was found with a significant separation from the hemisphere primary vertex, an 
artificial neural network was used to further separate b decays from charm and light quark background. 
This neural network has five inputs, derived from decay length, vertex multiplicity and invariant mass 
information. The vertex tag variable B for each hemisphere was then derived from the largest neural 
network output from any jet in the hemisphere |^. 

The distribution of the tagging variable B for 1994 data is shown in Figure |^(a-b) together with the 
expectation from Monte Carlo simulation. The distributions are shown separately for the barrel region 
(Icos^tI < 0.8) and the forward region (Icos^rl > 0.8) where the tagging performance is reduced due 
to the silicon microvertex detector acceptance. Three classes of hemisphere tags were defined: a 'tight' 
vertex tag T for hemispheres with B > 2.3, a 'medium' tag M for 1.6 < B < 2.3 and a 'soft' tag S for 
1.2 < S < 1.6. The hemisphere b-tagging efficiencies of the T, M and S tags are about 18%, 6% and 

5 %, and the fractions of tagged hemispheres originating from non-bb events are about 3 %, 15 % and 
25 %. The tagging efficiencies vary by up to about 15 % from year to year due to the differing silicon 
microvertex detector configurations. Events with B < have secondary vertices which are displaced 
from the primary vertex in the opposite direction to that of the jet momentum. The rate of these 
'backward tags' is sensitive to the detector resolution, and is slightly higher in data than in Monte 
Carlo, in both barrel and forward regions. The effect of this resolution mis-modelling is discussed in 



section 7.2 



Electrons and muons with momentum p > 2 GeV and transverse momentum pt > 1 GeV were also 
used to tag bb events. Electrons were identified in the polar angle region |cos0| < 0.96 using the 
neural network algorithm described in The identification relies on ionisation energy loss (dE/dx) 
measured in the tracking chamber, together with spatial and energy-momentum (E/p) matching 
between tracking and calorimetry. Photon conversions were rejected using another neural network 
algorithm Muons were identified in the same polar angle region by requiring a spatial match 
between a track reconstructed in the tracking detectors and a track segment reconstructed in the 
external muon chambers, as in |17|. 

The tagged lepton hemispheres were further enhanced in semileptonic b decays by using informa- 
tion from the lepton p and pt and its degree of isolation from the rest of the jet in a neural network 
algorithm Q|. The distribution of the neural network output variable {£ = e, /x) for identified 
electrons and muons is shown in Figure |l|(c) and (d). Hemispheres were defined to be tagged with the 
lepton tag L if any lepton in the hemisphere had an output Si > 0.7, corresponding to a b-tagging 
efficiency of 8% and a non-bb impurity of 20%. If a vertex T, M or S tag was also present in the 
hemisphere, it was ignored and the hemisphere considered only as an L tag. 

Events containing at least one hemisphere with a T, M, S or L b-tag were considered selected 
and used for the asymmetry analysis. Each combination of tags in the two hemispheres (T- nothing, 
L-nothing, T-T, T-S etc.) defined a separate tagging class, making a total of 14 tagging classes. In the 
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Figure 1: Distributions of (a-b) the vertex tagging variable B in different regions of 1008 6*^1; (c) the 
electron tag neural network Se', (d) the muon tag neural network for the 1994 data (points) and 
Monte Carlo simulation (histogram). The contributions from hemispheres containing light (uds), c 
and b quarks are indicated, and the cuts defining the hemisphere T, M, S and L tags are shown by 
the dashed lines. Distributions for the other years are similar. 
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data, 520133 b-tagged events were selected in one of the 14 classes, with a bb event tagging efficiency 
of about 54%. The most important tagging classes are those with a T, L, S or M tag opposite an 
untagged hemisphere, which comprise 33 %, 22 %, 14% and 13 % of the tagged data sample. 



5 Tagging the b production flavour 

The production flavour (b or b) of the b quark was determined independently in the two hemispheres 
of each selected event, irrespective of which hemispheres were tagged by the T, M, S or L b-tags 
described above. Up to four pieces of information were used per hemisphere: the momentum-weighted 
average track charge with two different weighting factors (evaluated for the highest energy jet in each 
hemisphere), the charge of a secondary vertex reconstructed in the hemisphere, and the charge of any 
kaon in the hemisphere, identified using dE/dx information. The jet charges can be calculated for every 
hemisphere, whilst the vertex and kaon charges are only available for a subset of hemispheres. The 
available information was combined using a neural network algorithm to produce a single production 
flavour tag variable Q for each hemisphere. Note that although semileptonic b decays were used to tag 
bb events, the charge of the lepton in L-tagged hemispheres was not used in the production flavour 
tag Q, in order to reduce the correlation with the b quark asymmetry analysis based on leptons Q. 
The jet charge Q-^^ was calculated for the highest energy jet in each hemisphere as: 



where p[ is the longitudinal momentum component with respect to the jet axis and qi the charge (±1) 
of track i, and the sum was taken over all tracks in the jet |19|. Two jet charges Q'^^^'^ and Q'^^^'^ 
were calculated for each jet, with the exponent k set to 0.5 and 1.0. The value k = 0.5 optimises the 
separation between hemispheres containing b and b quarks for a single jet charge @, and including 
the second jet charge with k = 1.0 provides a small amount of additional separation power, although 
the two jet charges are strongly correlated. 

For hemispheres containing a reconstructed secondary vertex, the charge of this vertex Qvtx was 
calculated as: 



and the uncertainty ag^ 



where Wi is the weight for track i to have come from the secondary, rather than the primary, vertex 
and the sum was again taken over all tracks in the jet. The weights Wi were obtained from a neural 
network algorithm using as input the track momentum, transverse momentum with respect to the 
jet axis, and impact parameters with respect to the reconstructed primary and secondary vertices, as 
in @. A well-reconstructed (small cTQ^t^) vertex charge close to -|-1 (—1) indicates a B"*" (B~) hadron, 
tagging the hemisphere as containing a b(b) quark, whilst a vertex charge close to zero indicates a 
neutral b hadron, {e.g. B" or B^), giving no information on the b quark production flavour. A vertex 
charge with large crg^tx cannot distinguish between charged and neutral b hadrons, and also provides 
no information on the b quark production flavour. 

Charged kaons produced from the b hadron decay can also be used to tag the b production flavour, 
via the underlying quark decay cascade b — > c ^ s. Candidate kaon tracks were selected in the highest 
energy jet in each hemisphere by using dE/dx information, requiring the track to have a probability to 
be consistent with a kaon of at least 5 %, and rejecting any tracks with a probability to be consistent 
with a pion exceeding 1 %. If more than one track in the jet was selected, the one with the highest 
weight Wi to come from the secondary vertex was retained. If no secondary vertex was reconstructed 
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in the jet, the weights were calculated using only the track momentum, transverse momentum and 
impact parameters with respect to the primary vertex. 

The hemispheres were then categorised into one of four classes, as follows: (1) neither vertex nor 
kaon charge, (2) vertex charge only, (3) kaon charge only (4) both vertex and kaon charges. The 
available tagging variables were combined using a neural network with up to five inputs: the two jet 
charges Q^^^'^ and Q^^^'^, the vertex charge Qvtx and error uq^^.^ and the weight Wi of the kaon 
track, signed by its charge. Separate neural networks were trained for each of the four classes. The 
continuous tagging variable Q is derived from the output x of the neural network, and is defined such 
that 

^ N^{x) + N^{x)' 

where A''b(x) and N^{x) are the number densities of Monte Carlo b and b hemispheres with a particular 
value of X. Hemispheres with Q = +1 (—1) are tagged with complete confidence as being produced 
from b(b) quarks, and hemispheres with Q = are equally likely to be from either. The modulus \Q\ 
satisfies \Q\ = 1 — 2^ where ^ is the 'mis-tag' probability, i.e. the probability to tag the production 
flavour incorrectly. The effects of B'^ and Bg mixing contribute to the mis-tag probability, since the 
decay flavour of a mixed b hadron is opposite to its production flavour. A similar tagging procedure 
was used in though with leptons rather than charged kaons included in the tagging information. 

The jet and vertex charge distributions are not charge symmetric, since detector effects cause 
differences in the rate and reconstruction of positive and negative tracks. These effects are caused by 
hadronic interactions in the detector material and the Lorentz angle in the tracking chambers [pi]] . 
They were removed by subtracting the mean value of each charge variable from the measured value 
before the calculation of Q. A final offset was then subtracted from Q as part of the asymmetry fit 
procedure. 

The distributions of Q in the four tagging classes are shown for all selected events in the 1994 data 
and Monte Carlo simulation in Figure ^. Some discrepancies are visible, particularly in class 1. These 
are not important since the tagging power of Q is measured directly from the data for bb events, 
and the corresponding effect on charm and light quark events is covered by the physics simulation 
systematic uncertainties. The four classes comprise about 25%, 36%, 15% and 24% of the data 
sample, and have effective mistag fractions of 33.1 %, 31.2 %, 32.5 % and 29.1 %0. 



6 Fit and Results 

The procedure which is used to derive the forward-backward asymmetry Ap-^ follows closely that in |Q , 
the main difference being the definition of the production flavour tag Q. In Q, Q was defined to be a 
jet charge with k = 0.5 and the vertex charge information was used in a separate fit. In this analysis, 
Q is defined to be the output of the artificial neural network tag described in the preceding section, 
which incorporates all information from the jet, vertex and kaon charges. 

In the case of a data sample consisting of only bb events without contamination from lighter 
quarks, and neglecting acceptance effects, it can be shown that 

{Qf - Qb) = ^FB • '^b , 

where Qp and Qb are the production fiavour tags in the forward and backward hemispheres as mea- 
sured in the data, and {Qp — Qb) is the mean difference of the production flavour tag values in the two 
hemispheres. The variable (5b is equal to {Q^ — Q+), where Q- (<5+) is the charge in the hemisphere 
containing the negatively (positively) charged primary b (b) quark. It measures the mean charge sep- 
aration between negatively and positively charged hemispheres and is sensitive to the tagging power 



of g m 



■^The effective mistag fraction measures the fraction of hemispheres that are incorrectly tagged, after weighting to 
take into account the confidence with which a hemisphere is tagged as b or b, measured by the value of \Q\. 
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Figure 2: Distributions of the production flavour tag Q in the four tagging classes for 1994 data 
(points) and Monte Carlo simulation (histogram). The contributions from hemispheres containing b 
and b quarks are shown, together with the charm and light quark (uds) background. Distributions 
for the other years are similar. 
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In the presence of charm and hght quark backgrounds in the data sample, and in case of tagging 
efficiencies varying as a function of | cos 9t\, the above equation must be modified: 

(Qf-Qb)= E ^fPfCfSf4^. (1) 

flavours / 

In this equation, Sf is +1 for down-hke and —1 for up-hke quarks, Pj is the fraction of events of 
flavour / in the data sample, derived mainly from the data as described below, and Sf is the charge 
separation for flavour / which is determined directly from the data for bb events and from Monte 
Carlo simulation for the other quark flavours. The factors Cf, which are taken from Monte Carlo 
simulation, account for variations of the asymmetry and the tagging efficiency with cos 9t and are 
given by 

_8 Jr]f{y)ydy _8 Eevcnts 1/ .o^ 
^ 3jrjf{y){l+y^)dy 3 Eevents(l + ' ^' 

with y = I cos^tI; Ty^ is the efficiency to tag an event of flavour / in a given small interval of y. The 
sums run over all tagged Monte Carlo events of flavour /. 

In the absence of correlations between the charges Q measured in both hemispheres (except for 
primary quark charges), and if there are no charge biases, i.e. if the mean hemisphere charge of all 
hemispheres is zero, the mean charge separation 6 can be derived directly from the data according to 

j = -{Qf-Qb). (3) 
In the presence of charge biases (due to detector effects as discussed in Section Isl) and correlations 



between the hemispheres, equation (^) must be modified and becomes |22| 



6^ _ -(QF-QB)+/0-g-^ + M^ 
4 ~ l + p 

where fj, and are the mean and variance of the hemisphere charge for all hemispheres measured 
from data. The variable p is the correlation between Q_ and Q^, evaluated in Monte Carlo bb events 
and given by 

{Q+ ■ Q.) - • 



1-2 f-iT-i/^ ^2 



(4) 



where ctq^ and <Tq_ are the variances of the distributions of Q+ and Q-. In the presence of charm 
and light flavour background, the measured charge separation 6 receives contributions from all flavours 
according to their fraction of the data sample: 

flavours / 

The values of 6f for charm and light quarks can be derived from Monte Carlo simulation according to 

Sf = Sf (Q_ - Q+)f 

where {Q- — Q+)f is determined in events of flavour /. This allows 5b to be determined from 
equation once the flavour fractions Pf are known, and then allows the b quark asymmetry A^-q to 
be derived by solving equation (||), using assumed values for the charm and light quark asymmetries. 

These asymmetries, together with fractions Rf of hadronic Z" decays to each quark flavour, were 
set to the Standard Model expectations as calculated by ZFITTER 6.36 |23|, and are given in Table |. 



Since part of the data sample was taken at centre-of-mass energies above and below the TP peak, the 
values of {Qf — Qb) were calculated separately for three energy bins, allowing the corresponding values 
of Apg to be determined. The energy bin limits and mean centre-of-mass energy in each bin are also 
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Flavour 


Rf 


4b ((\^) =89.50 GeV) 


4b ((\^) =91.26 GeV) 


4b ((\^) =92.91 GeV) 






(88.40-90.40 GeV) 


(91.05-91.50 GeV) 


(91.70-94.00 GeV) 


bb 


0.2155 








cc 


0.1726 


-0.0309 


0.0633 


0.1211 


ss 


0.2196 


0.0595 


0.0964 


0.1189 


uu 


0.1728 


-0.0308 


0.0632 


0.1207 


dd 


0.2196 


0.0595 


0.0964 


0.1189 



Table 1: Values for Rf for bottom, charm and light quarks and forward-backward asymmetries for 
charm and light quarks as calculated from ZFITTER [^] for the three different energy bins. For each 
energy bin, the mean energy and the bin limits are also given. 



ej' is the efficiency for an event of flavour / to be tagged by tag k in one hemisphere and tag I in 
the other hemisphere. Deviations of D¥ from unity account for the fact that the tagging in the two 



given in Table |^. All data were used for the calculation of 5^ since it does not vary significantly with 
centre-of-mass energy. 

The above procedure was applied separately to each of the 14 tag classes and in five bins of |cos 9t\ 
with bin edges at |cos0t| = 0, 0.2, 0.4, 0.6, 0.8 and 0.95. The flavour fractions Pf were determined 
simultaneously for all 14 tag classes in each Icos^^l bin using the measured tagging rates for each 
class. Within the tagging class kl, where k,l={T, M, S, L or nothing }, i.e. tagged in one hemisphere 
by b-tag k and in the other hemisphere by b-tag /, the fraction of events Pj' of flavour /={b,c,uds} 
is given by 

where the sum in the denominator runs over b, c and light quark (uds) flavours, is the hemisphere 
tagging efficiency of tag k for flavour /, and the correlation is defined by = ej'/(£/£/)) where 

hemispheres is not completely independent. For light quark events, these correlations have negligible 
effect and are set to one. The fraction of hemispheres fl in the |cos 9t\ bin that are tagged by b-tag i, 
and the fraction of events /J^' tagged by b-tag k in one hemisphere and b-tag / in the other hemisphere 
are then given by 

fl = eL-^b + Cc-^c + euds(l ~ -^b - -Rc) , 
fd = ^h^h^h Rh + e^e^D^ Rc + euds^^uds(l - -^b - Rc) ■ 

This system of equations was solved using a maximum likelihood fit. The uds tagging efficiencies 
and all correlation terms were fixed to values determined from Monte Carlo simulation, and the b- 
and c-tagging efficiencies of the T, M, S and L tags were allowed to vary in order to minimise the 
difference between the observed and predicted tag fractions. Then the flavour fractions Pj' for each 
tagging class were calculated from equation defining the untagged efficiencies for each fiavour as 
^nothing _ _ gT _ gM _ gS _ gL_ 'pj^g values for R\^ and Rc were computed using ZFITTER (see 
Table!). 

The asymmetry fit procedure was applied separately to the events collected in each of the five 
I cos I bins and 14 tag classes, resulting in 70 different values for Ap-^ for each energy bin and 
data taking period (1991-2, 1993, 1994, 1995 and 1996-2000). All the measurements for each 
energy point were then combined, weighted according to their statistical errors. The mean charge fiow 
{Qp — Qb) as a function of |cos 9t\ is shown for peak events with either hemisphere tagged by a L, 
S, M or T tag in Figure |3|(a-d), and for all tagged events at the off-peak energy points in Figure ^(e) 
and (f). The expected distribution from the result of the asymmetry fit is also shown in each case. 
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Figure 3: The mean charge flow as a function of |cos0t| in the entire dataset (points with error 
bars showing the statistical errors) , together with the prediction from the asymmetry fit (histogram) . 
Figures (a-d) show the distributions for Z° peak events tagged in either hemisphere with an L, S, M 
or T tag, whilst (e) and (f) show all tagged events at the two off-peak energy points. 
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These fitted asymmetry values do not correspond directly to the b quark forward-backward asym- 
metry because of the effects of gluon radiation from the primary quark pair and the approximation 



of the original quark direction by the experimentally measured thrust axis |24|. The effects of gluon 
radiation have been calculated to second order in as, using the parton level thrust axis to define the 
asymmetry [^5|. The correction needed to go from the parton level to the hadron level thrust axis 
(calculated using all final state particles without detector effects) has been determined using Monte 
Carlo hadronisation models However, these corrections cannot be applied directly to this anal- 
ysis, since the determination of the production flavour tagging power directly from the data already 
accounts for most of the effects of gluon radiation. Therefore, a large sample of Monte Carlo simulated 
events was used to determine the correction to the measured asymmetries directly, by comparing the 
asymmetry fit results on this sample with the true primary b quark asymmetry. This correction factor 
was rescaled so as to take the quark to hadron level correction from the theoretical calculation dis- 
cussed above, since the calculation is expected to be more accurate than the Monte Carlo simulation 
for this part of the overall correction. Combining all effects, the measured asymmetries were scaled 
by 0.9923 it 0.0063 it 0.0038 to determine the quark level asymmetries, where the first error is due 



to theoretical uncertainties in the quark to hadron level calculation |24,26| and the second to Monte 
Carlo statistics. 

The measured asymmetry values, after correcting to the quark level, are: 

A^^ = 0.0582 ± 0.0153 at yi = 89.50 GeV 
A^B = 0.0977 ± 0.0036 at ^ = 91.26 GeV 
A^B = 0.1221 ±0.0123 at ^ = 92.91 GeV 

where the errors are statistical only. 



7 Systematic errors 

Systematic errors on Ap-Q arise from uncertainties in the input quantities which are taken from Monte 
Carlo, namely the b-tagging correlations Dy^, the production flavour tag correlations /Jb, the charge 
separations df for charm and light quark events, the light quark tagging efficiencies euds and the 
efficiency correction factors Cj. Additional uncertainties result from material asymmetries in the 
detector and the calculation of the QCD corrections to the raw measured asymmetries. The systematic 
errors are summarised in Table § and discussed in more detail below. 

7.1 Hemisphere correlations 

The following uncertainties in the simulation of b hadron production and decay affect the estimates of 
both the b-tag correlations Db ™d the production flavour tag correlations pb- They were estimated 
by reweighting the Monte Carlo sample used to derive the correlations and repeating the asymmetry 
fit using the modified parameters. 

b quark fragmentation: The Monte Carlo was reweighted so as to vary the average scaled energy 
{xe) of weakly decaying b hadrons in the range {xe) = 0.702ib0.008, as determined by the LEP 



electroweak working group |26|. The fragmentation functions of Peterson et al., Collins and 
Spiller, Kartvelishvili et al. and the Lund group [|l^,^ were each used as models to determine 
the event weights, and the largest observed variations in Apg were assigned as the systematic 
errors. 

b hadron production fractions: The fractions of b quarks hadronising to form Bg mesons and Ab 
baryons were varied in the ranges /(b B^) = (10.7±1.4) % and /(b Ab) = (11.6±2.0) % [||. 

b lifetimes: The lifetimes of b mesons were varied by ±0.02 ps and b baryons by ±0.05 ps, based on 
the uncertainties on the measured values |28]. 
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0.00044 


0.00047 


0.00082 


Silicon alignment 
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0.00003 


0.00009 


Electron fake rate 


0.00002 


0.00014 


0.00021 


Muon fake rate 


0.00004 


0.00022 


0.00034 


Material asymmetry 


0.00009 


0.00009 


0.00009 


Total detector 


0.00087 


0.00093 


0.00117 


QCD and thrust axis correction 


0.00043 


0.00073 


0.00091 


Event selection bias 


0.00002 


0.00011 


0.00016 
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Table 2: Summary of systematic uncertainties on the measured values of Ap^ for the three different 
energy bins. 
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b decay charged multiplicity: The average charged decay multiplicity of b hadrons was varied by 
±0.062, reflecting the accuracy of the measurements by LEP experiments 

Hemisphere b-tagging correlations^ were extensively studied for the measurement of i?b P and 
found to result from two main sources: (i) kinematical correlations between the momenta of the two b 
quarks in a bb event (due to hard gluon radiation and soft particles produced in the hadronisation of 
the two b quarks); and (ii) geometrical correlations caused by the strong dependence of the b-tagging 
efficiency on | cos 0^1- Kinematical correlations were found to be small, and modelled in the Monte 
Carlo to a precision of = 0.0022 [Q. The systematic error for the analysis was assessed by 
simultaneously changing the b-tag correlations L'b for all 70 analysis bins by ±0.0022 and repeating 
the asymmetry fit. 

The geometrical correlation contributions to Dy, were measured directly from the data in each of 
the 70 analysis bins, using the rate of tagged events as a function of the polar and azimuthal angles 
of the thrust axis, as described in |5|. At high |cos0r|) the geometrical correlations are large, rising 
to around Db — 1 = ±0.5 for the T-T double-tagged events at |cos^r| > 0.9. Here, the assumption 
of independent geometrical and kinematic correlations breaks down, and their sum overestimates the 
overall correlation determined directly from the Monte Carlo single and double tag efficiencies. To 
account for possible Monte Carlo mis-modelling of this effect, the full difference between the correlation 
component sum and the overall correlation was taken as an additional systematic error on D^,, for all 
analysis bins where this difference was statistically significant. To assess the corresponding systematic 
error on Ap^, the asymmetry fit was repeated with the D\, correlations for all such bins shifted 
simultaneously. 

The largest single contribution to the overall systematic uncertainty arises from the correlation 
Pb between the production flavour tag determinations in the hemispheres of the primary b and b 
quarks, Q- and (see equation (^)). The values of pb are typically between zero and —10%, 
depending on tagging class and |cos 6t\- The origin of the correlation is primarily events with significant 
gluon radiation, which reduces the momenta of both b hadrons, and may also lead to a third jet 
shared between the two hemispheres. Both of these effects reduce the average tagging power in both 
hemispheres of the event, leading to a tagging correlation. This can be seen in Figure |^(a), which shows 
the overall correlation /Ob averaged over all 70 tagging bins as a function of thrust. The correlation is 
stronger in events with low thrust values, corresponding to significant gluon radiation and a three-jet 
topology. 

In order to check the modelling of the correlation pb in the Monte Carlo simulation, the variables 
that contribute to the correlation were determined, and studied using a reweighting procedure. For a 
variable X possibly contributing to pb, the correlation 

was determined in Monte Carlo events, where X+ and X- are the values of the variable in the 
hemispheres of the primary b and b quarks. Event weights were then determined in small bins of 
and X- such that the correlation between and X_ could be removed, and the effect on the 
overall correlation pb was studied. Using this procedure, the two jet charges Q^^^'^ and Qf^^'^ and 
the vertex charge Q^tx were identified to be the relevant variables contributing to the non-zero value 
of Pb- 

Since the determination of Cx is possible only in Monte Carlo events, where the hemisphere of 
the original b quark is known, an additional set of correlations were defined according to 

~ {\Xp\-\Xb\)-{\Xp\)-{\Xb\) 

L-X = , 

C^XpO-Xg 

^The hemisphere b-tagging correlations Dh were denoted by in reference fej. 
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Figure 4: Studies of the production flavour tag liemisphere correlations pb (see text): (a) Average 
value of pb ^ ^ function of thrust; (b-d) correlation test variables C for Q'^^^-^, Qf^^'^ and Qvtx in 
data (with shaded error band) and Monte Carlo (as a function of correlation reweighting factor a), 
together with corresponding values of pbi (e) correlation between Q^tx and production flavour tag Q in 
data and Monte Carlo (as a function of correlation reweighting factor a). The lines through the points 
in Figures (b~e) show the result of linear fits used to parameterise the dependence of the correlations 
on a. 
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where Xp and Xb are the values of the variable X in the forward and backward hemispheres. These 
Cx are sensitive to the correlation Cx and can also be measured in data, to check the Monte Carlo 
modelling of the hemisphere correlation of the variable X. This was done by applying scale factors a 
to the weights used to remove the correlation, and calculating Cx in Monte Carlo events as a function 
of a. The values of Cx were then compared to Cx measured in data to determine the range of a 
allowed by the data. This procedure is illustrated in Figure Q(b-d), which show for the three variables 
of interest the data value of Cx together with the Monte Carlo values of Cx and the overall correlation 
Pb as a function of the scale factor a. 

For the two jet charges Q^^^'^ and Q^^"^'^ ■, the agreement between data and unweighted Monte 
Carlo is good, and the correlations Cx are sensitive measurements of the corresponding contributions 
to /9b- The systematic error on A^-q was calculated by varying the amount of reweighting in the 
range allowed by the data, recalculating the values of pb for all 70 analysis bins and refitting the 
asymmetries. This procedure results in systematic uncertainties of 0.00048 and 0.00070 on Apg in the 
Z peak energy bin for the jet charges Qjet and Qjet respectively. For the vertex charge (c^vtx? 
Figure ^d) shows that data and Monte Carlo are consistent (within two standard deviations), but that 
the size of the statistical error is too large and the effect on pb too weak to draw precise conclusions 
about the Monte Carlo modelling of the Q^tx correlation. Figure ^(e) shows a more sensitive variable, 
the correlation C of the vertex charge in one hemisphere and the production flavour tag in the other 
hemisphere, where again good agreement between data and Monte Carlo is seen, giving confidence 
in the Monte Carlo description of the vertex charge correlations. However, since this test is more 
indirect and cannot be directly applied to the vertex charge reweighting, the systematic error on Apg 
is assessed by reweighting based directly on the vertex charges in the two hemispheres, so as to increase 
or decrease the correlation Cq^^^ by 50 %. This gives an additional systematic error on of 0.00026. 
Adding the contributions of the three variables in quadrature gives a total error due to the modelling 
of the pb hemisphere correlations of 0.00089. 

7.2 Detector Simulation 

Both the tagging correlations and pb, and the other input parameters 5, Cuds and Cf are sensitive 
to details of the detector simulation, in particular the tracking and lepton identification performance. 

Tracking resolution: The error due to uncertainties in the tracking resolution was assessed by ap- 
plying a global 10 % degradation to the resolution of all tracks, independently in the r-cp and 
r-z planes, as in This degradation accounts for the discrepancies between data and Monte 
Carlo backward tagging rates seen in Figure ^(a) and (b). The resolution is also sensitive to the 
efficiency for associating silicon hits to tracks, which was varied by ±1 % in the r-(j) and ±2% 
in the r-z planes to cover residual discrepancies between data and Monte Carlo hit association 
rates. 

Silicon alignment: The b-tagging performance is sensitive to knowledge of the radial positions of 
the silicon microvertex detector wafers, which are known to a precision of ±20 pm from studies 
of cosmic ray events The corresponding uncertainty is assessed by displacing one or both 
barrels radially by 20 pm in the simulation. 

Lepton identification: The light quark tagging efficiencies e^ds for the L tag are sensitive to the 
number of charged hadrons mis-identified as electrons and muons, which is modelled in the 
Monte Carlo to precisions of ±21 % for fake electrons and ±9 % for fake muons The analysis 
is not sensitive to the Monte Carlo description of the efficiencies for identifying real leptons, 
since these occur primarily in bb and cc events, whose tagging efficiencies are measured directly 
from the data. 
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7.3 Charm and light quark physics 

The following uncertainties in the simulation of charm physics affect the asymmetry analysis through 
the input values of 5^, the charge separations in cc events. The effect of corresponding uncertainties 
on the charm hemisphere b-tagging and production flavour tagging correlations is negligible. 

Charm fragmentation: The Monte Carlo was reweighted so as to vary the mean scaled energy 
{xe) of charm hadrons in cc events in the range {xe) = 0.484 it 0.008 [26|, using the same four 
fragmentation functions as for bb events. 

Charm hadron production fractions: The production fractions of the weakly decaying charm 
hadrons were varied according to the measurements performed at LEP ||2^, as averaged by 
the LEP electroweak working group The contribution from was scaled by 1.15 it 0.05 
to account for other weakly decaying charm baryons. The dependence on the production rates 
of excited charm states was determined by varying the fractions of charm quarks hadronising to 
produce D*+, 0*° and D* in the ranges 0.239 ±0.007 0.218 ±0.071 [||] and 0.13 ±0.13 H] 
whilst keeping the production fractions of weakly decaying charm hadrons constant. The depen- 
dence on the production of orbitally-excited charm states (D**) was found to be negligible. 

Charm hadron lifetimes: The lifetimes of the weakly decaying charm hadrons were varied sepa- 
rately according to the measured values [p^. 

Charm hadron decay multiplicities: The average charged hadron decay multiplicities of D+, D*^ 
and mesons were varied according to the measurements of MARK III [^]. The charged 
decay multiplicity of charm baryons, for which no measurements are available, was varied by 
±0.5. The number of vr*^ produced in D meson decays were also varied according to the available 



measurements |32|. The branching ratio of charm hadrons to long-lived neutral strange particles 



(K*^ and A) were varied according to the uncertainties quoted in |28|. In each case, the other 
branching ratios were held constant whilst the branching ratio in question was varied. 

Charged kaon production in charm decays: The tagging performance of the production flavour 
tag in cc events is sensitive to the number of charged kaons produced in D meson decays. These 



were varied according to the measured values |28| 



Uncertainties in the simulation of light quark events affect both the tagging efficiencies euds a-nd 
the charge separations (5u, 5 a and 5^^. The inclusive production rates of K'' mesons and A and other 
weakly decaying hyperons were varied in the Monte Carlo by ±3.4 %, ±6.5 % and ±11.5 % respectively. 



corresponding to the precision of the OPAL measurements [33| combined with an additional uncer- 



tainty to take into account the extrapolation of the inclusive production rates to light quark events 



only. Additionally, HERWIG 6.2 |3j and ARIADNE 4.08 |§5| were used as alternative fragmentation 
models for the simulation of light quark events. 

The production of heavy quarks via the gluon splitting processes g — > cc and g ^ bb affects 
the properties of both charm and light quark events. The rates were varied independently in the 
ranges /(g ^ cc) = (2.96 ± 0.38) % and /(g ^ bb) = (0.254 ± 0.051) % according to LEP and SLD 
measurements l2f 



7.4 Other uncertainties 

The parameters Cj which correct for the variation of tagging efficiency with Icos^j-I within each 
bin are calculated separately for each flavour using Monte Carlo simulation (see equation The 
simulation was checked by studying the rate of tagged events as a function of |cos^t| within each of 
the 70 analysis bins, and reweighting the Monte Carlo b-tagging efficiency in small bins of |cos0t| to 
reproduce the data distribution. The flavour dependence of the efficiency corrections was checked by 
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Derivative 


(^) = 89.50 GeV 


{^/s) = 91.26 GeV 


(V^) = 92.91 GeV 




-0.0175 


-0.0166 


-0.0168 


d^pg/ d^pg 


0.0221 


0.0215 


0.0222 


d^Fs/ d^FB 


-0.0222 


-0.0224 


-0.0222 


d^^^/dA^B 


0.0686 


0.0695 


0.0694 




0.0032 


0.0215 


0.0331 




0.0000 


0.0000 


0.0000 


dA%/dRc 


0.0030 


0.0134 


0.0202 


dA^B/di?b 


-0.0692 


-0.3482 


-0.5104 



Table 3: Derivatives of ^pB with respect to the assumed Standard Model parameters, for each energy 
bin. The value of Rd is constrained to 1 — i?b — -Rc — -Ru — -Rs- 



setting all charm and light quark Cf parameters to the corresponding Cb. The resulting changes in 
the fitted values of A^^ were negligible in both cases. 

As discussed in section the jet and vertex charge distributions are not charge-symmetric due 
to detector effects. A difference in the amount of detector material in the forward and backward 
hemispheres could lead to different offsets in the two hemispheres and bias the measured value of ^pB ■ 
This material asymmetry was measured by studying the rate of identified photon conversions as a 
function of cosOt- For |cos0t| < 0.8, the conversion asymmetry was found to be consistent with zero 
to a precision of ±0.3 %. For |cos 9t\ > 0.8, the backward hemisphere was found to contain 7.2 ±0.5 % 
more material than the forward hemisphere, due to readout electronics and cabling. The systematic 
error on Ap-g was calculated by assuming that all of the jet charge offsets are due to material effects 
and differentially shifting them in the forward and backward hemispheres according to the measured 
material asymmetries as a function of |cos0t|- 

The correction from the asymmetry measured using the experimental thrust axis to the primary 
b quark asymmetry is known to a precision of 0.74%, including both theoretical uncertainties and 
Monte Carlo statistics. Note that the full size of the theoretical error is used, even though part of the 
gluon radiation correction is absorbed by the determination of the tagging power from the data, as 
discussed in Section ^. 

The hadronic event selection requirements are (0.25 ± 0.15) % more efficient for bb events than 
for charm and light quark events ||5|, due mainly to the requirement of at least seven charged tracks. 
This leads to a small uncertainty on the flavour composition of the sample, and a corresponding error 
on ^pB- The LEP centre of mass energy is known to a precision of 18MeV for the peak running 
in 1992, around 5MeV for 1993-95 and 12MeV for the Z° calibration data taken in 1996-2000 
Taking year-to-year correlations into account, and assuming the Standard Model dependence of Ap-^ 
on y/s, this leads to the uncertainties given in Table § on the asymmetries at the quoted values of ^/s. 

The asymmetries measured in each year of the data, and in different tag classes and |cos0t| bins 
are consistent. The results were found to be stable when the b-tagging cuts defining the different 
tag classes T, M, S and L were varied, and no additional systematic error is assigned. The fit was 
also tested on large samples of simulated Monte Carlo events with different true asymmetry values, 
and no evidence of a bias was seen. As a further cross-check, the rates of single- and double-tagged 
events were used to measure i?b as a function of |cos 9t\ and tag type, and results consistent with the 
world-average value of R^ were obtained in all cases. 

The fitted asymmetry values depend on the assumed values for the fraction of hadronic decays 
to each quark flavour, and the charm and light quark asymmetries at each energy point. The values 
used are given in Table and the derivatives of the measured asymmetries with respect to these 
parameters are given in Table |^ Taking the values of i?b -Rc and Ap-^ from measurements |28| rather 
than ZFITTER would result in additional uncertainties on ^^b at = 91.26 GeV of 0.00026, 0.00006 



20 




Figure 5: The measured b quark asymmetry as a function of energy (points with error bars showing 
both statistical and systematic errors), together with the Standard Model expectation calculated using 
ZFITTER m, with the Higgs mass varied between 114 GeV and 1000 GeV. 



and 0.00031 for i?b, Rc and Ap-^ respectively. 



8 Conclusions 

The b quark forward-backward asymmetry has been measured using an inclusive tag at three energy 
points around the peak. The results, corrected to the primary quark level, are: 

A^^ = 0.0582 ±0.0153 ±0.0012 at ^ = 89.50 GeV 
A^B = 0.0977 ± 0.0036 ± 0.0018 at ^ = 91.26 GeV 
yl^B = 0.1221 ±0.0123 ±0.0025 at ^ = 92.91 GeV 

where in each case the first error is statistical and the second systematic. The results are shown as a 
function of ^/s in Figure ^, together with the Standard Model expectation calculated using ZFITTER 
[ p^ . Using the ZFITTER prediction for the dependence of A^-q on ^/s, the three measurements 
are shifted to mz (91.19 GeV), averaged and corrected for initial state radiation, 7 exchange, j — 
interference and b quark mass effects. The resulting value for the Z^ pole asymmetry Ap^ is 

Ap'° = 0.1002 ± 0.0034 ± 0.0018 

where again the first error is statistical and the second systematic. Within the framework of the 
Standard Model, this corresponds to an effective weak mixing angle for electrons of 

sin^ e^'"" = 0.23205 ± 0.00068 . 

This result is one of the most precise measurements of the b quark forward-backward asymmetry to 
date. It is in agreement with, and supersedes, the previous OPAL result using jet and vertex charge |^, 
and is also in agreement with the OPAL measurement using leptons and other measurements of 
A^^ at LEP g. 
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